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Abstract 



In this paper, we investigate the generahzed Saez-Ballester scalar-tensor theory of grav- 
£^ \ ity via Noether gauge symmetry (NGS) in the background of Bianchi type I cosmological 

i-J^ . spacetime. We start with the Lagrangian of our model and calculate its gauge symmetries 

and corresponding invariant quantities. We obtain the potential function for the scalar field 
in the exponential form. For all the symmetries obtained, we determine the gauge fun- 
tions corresponding to each gauge symmmetry which include constant and dynamic gauge. 
We discuss cosmological implications of our model and show that it is compatible with the 
observational data. 
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I. INTRODUCTION 

Several cosmological observations indicate that the observable universe is undergoing a phase 
of accelerated expansion [l]. There are two major approaches to address the problem of cosmic 
acceleration: either introducing a 'dark energy' component in the universe and study its dynamics 

or interpreting it as a failure of general relativity (GR) and consider modifying GR theory, 
termed as the 'modified gravity' approach [3]. Although both approaches have novel features with 
some deep theoretical problems, we here focus only on the modified gravity approach. 

One of the earlier modifications to Einstein's general relativity was termed Brans-Dicke gravity, 
in which besides a gravitational part, a dynamical scalar field was introduced to account for a 
variable gravitational constant [4]. This modification was introduced due to lack of compatibility 
of Einstein's theory with the Mach's principle. Another cosmologically viable model is the covariant 
scalar-tensor- vector theory which allows the gravitational constant G [5]. The theory can explain 
successfully galaxy rotation curves and cluster data without non-baryonic dark matter. Moreover 
the theory is consistent with solar system observational tests. Later Saez & Ballester Q] introduced 
a scalar-tensor theory of gravity in which metric is coupled to a scalar field. Here the strength of 
the coupling between gravity and the field was governed by a parameter w. With this modification, 
they were able to solve a 'missing- mass problem'. Several aspects of Saez & Ballester theory in 
relation to Bianchi cosmological models have been explored in literature Q]- 

Noether symmetries are the symmetries of the Lagrangian. In literature, the approach of 
Noether symmetry is used to obtain exact forms of gravitational theories including f{T) gravity. 



where T is torsion scalar 
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^, f{R) gravity, R being Ricci scalar and scalar-tensor theories 



1C|. Spherically symmetric solutions in f{R) gravity via Noether symmetry were discussed in 



. This approach gives a power-law evolutionary form of scale factor which is consistent with 
the astrophysical observations. Moreover the evolution of state-parameter obtained in such an 
approach also gives a phantom crossing behavior of dark energy 12]. The Noether symmetry 



1^ 



approach has been applied to Bianchi cosmological models in literature: Capozziello et al 
investigated the Bianchi universes via Noether symmetries. Camci &: Kucukakca \\\ studied 
the Noether symmetries of Bianchi type I, III and Kantowski-Sachs spacetimes in scalar coupled 
theories. They obtained the exact solutions for potential functions, scalar field and the scale 
factors, see also [l^ which is a similar work as [14]. Scalar-tensor theories have been investigated via 
Noether symmetry [16] but via NGS approach, the analysis is not reported earlier. The application 
of Noether theorem in higher order theory of gravity turned out to be a powerful tool to find the 
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solution of the field equations and physically reasonable solutions like power law inflation have 
been discussed in The NS approach has been applied to pure gravity model with variable 



cosmological constant A and gravitational constant G in [1 

In this paper, we consider a Bianchi type I spacetime in the framework of Saez-Ballester theory 
of gravity. We set up a Lagrangian in which the metric variables and scalar potential play the role 
of dynamical variables. The Lagrangian is so constructed that its variation with respect to the 
metric components and the scalar potential leads to the correct equations of motion. We explicitly 
calculate the form of scalar potential by demanding the Lagrangian admits the desired Noether 
symmetry. Unlike the usual approach to obtain Noether symmetries CxL = as followed in ^, 
we employ the fuh Noether Gauge Symmetry Condition (X^^IL + (Dt£,)L = DtG{t,(l),A,B,C)) 
introduced earlier in [19]. The advantage of this later scheme is that it yields extra symmetries 
then the former one, hence full depth of the theory is realized in this manner. 

The plan of the paper is as follows: In section II, we model our system by writing the Lagrangian 
and then deriving the dynamical equations of motion for Bianchi type I spacetime. In section III, 
we consider pure vacuum solution and construct a Lagrangian. Using it, we solve system of coupled 
differential equations to obtain Noether gauge symmetries and corresponding invariant quantities. 
We discuss some cosmological implications of our model in section IV. We conclude in section V. 



II. THE MODEL 

The metric of Bianchi model of type I in coordinates x'^ = {t, x, y, z) is [2C 

g^, = diag{l, -A\t), -B\t), -C\t)) (1) 

The exact solutions of Einstein field equations based on metric ([T]) have been investigated in detail 
, 2l| . The geometrical quantities of ([l]) are the average scale factor a = ^/ ABC ; 



in the literature 
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the volume of the spacelike hypersurface, defined hy V = ABC = a?; the generalized (or mean) 
Hubble parameter: H = ^Yli^i where Hi = dtlog{Ai), Ai = {A,B,C}. For the isotropic case, 
A = B = C, the mean Hubble parameter converts to the Friedmann-Robertson- Walker form 
H = dt log(a) where a is the scale factor. In this paper we are interesting to investigate the 
anisotropic models in which the cosmology described by metric ([1]) with A ^ B ^ C . 

We consider the case of a homogeneous but anisotropic Bianchi type-I model with a scalar field 
4> based on a non-standard scalar-tensor theory. The action of this model reads [g] 

S = ^J V^d^x^R + |</<V,M</''^ - V {<!>)) , (2) 
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where k and oj are arbitrary dimensionless constants. Choosing k = reduces our model to the 
minimally coupled massless scalar field coupled to Einstein gravity. Different aspects of this model 
have been explored in the literature 



22[ |. Varying ([2]) w.r.t the metric g^u lead to a generalized 



Einstein equation 



The generalized Klein-Gordon equation for scalar field is 



Using metric ([I]) in the field equations ([3]) and (jl]), we obtain 
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E(4^+4^+4^4^) = V>', (5) 



^ ^ 



Ef^ = -^-0'=0^ + m. (6) 



The Klein-Gordon equation is 



III. NGS ANALYSIS 

We eliminate the terms Ai and obtain the following Lagrangian which is suitable for calculating 
the gauge symmetries: 

L{t,(l),A,B,Cj,A,B,C) = (^'^4>^^'^ -V{cl))^ ABC -2{ABC + ABC + AbC), u ^ (8) 

Varying the Lagrangian ([8]) w.r.t. (p, A, B and C, we get a system of Euler-Lagrange equations 
(or field equations): 

.. _ 2(I)ABCV^ + 2uj(j)^<j){kABC(j) + 0(iSC + ABC + ABC)) 



lujcpf'+^ABC 



(9) 



.. _ AjABC - ABC - ABC) + 2ABCV - w(j)''^'^ABC 

•• A{ABC - ABC - ABC) + 2ABCV - w(p''<j)^ABC ^ , 

^ = 8AC ' 

^ _ 4(iSC - iSC - ABC) + 2A5Cy - wcl)''(j)'^ABC 

^ - 8AB • ^^^^ 
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The Noether symmetry is given by 

where the coefficients ^, r/j, {i = 1,2,3,4) are determined from the Noether symmetry conditions. 
The first order prolongation of the above symmetry to the first-order jet space comprising of all 
derivatives is 

' dT dA dB dc 

The Noether gauge symmetry condition is jiol 

+ [Dt^L = DtG{t, 4>, A, B, C), (14) 

where G is the gauge function. We emphasize here the difference between the Noether and Noether 
gauge symmetries: In fact the NS is a very special case of NGS i.e. ignoring the gauge function 
and first prolongation, we find the restricted (or a special form of) Noether symmetry. The set of 
Noether symmetries is always a subset of Noether gauge symmetries. 

The condition (jl4p yields the following system of linear partial differential equations 

^0 = 0, 

Cr/3,A + Brji^A = 0, 
Cr]2,B + Arji^B = 0, 
Br]2^C + Ar]3^c = 0, 
Ga + 2Br]4,t + 2Crj3,t + VABCU = 0, 



Gb + 2Cr]2,t + 2Ar]4,t + VABC^b = 0, 
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Gc + 2Ar]3,t + 2Br]2,t + VABC^c = 0, 
2Cm,^ + 25?74,^ - ViA^^^^BG = 0, 
2Gr]2,^ + 2Ar]4^^ - r^i^BWcj)^ ABC = 0, 
2Bri2,4> + 2Ar,3^^ - rji^cw4>^ ABC = 0, 
A-ni,A + Bri4,B + C^t + C {r]3,B - ^t) + C (??2,A - 6) + ??4 = 0, 
Bm,B + Gr)2,c + ACt + A {ris^B - ^t) + A (r?4,c - Ct) +V2 = 0, 
Am,A + Grj3^c + Bit + B {7]2,A -it)+B (774,0 - 6) + ??3 = 0, 
A:^BC7r?i + </.(r/4A5 + r/s^C + 7?3^C7 + Ci^BC + 2 - 6) wABC) = 0, 
+ r]iABCV^ + V{ABCit + mAG + rj^AB + r/a^C) = 0, 



corresponding to the gauge functions G{t, </>, A, B, C). We numerically solve above system of linear 
partial differential equations. The potential function V{(j)) is an arbitrary function whose form will 
be specified by the determining equations. We have the following cases. 



The above system (|15p altogether yields nine Noether symmetries comprising of translation, 
scalings and other symmetries. 



G^ + ABC{V(^-m,tw^'') = 0. 



(15) 



A. T/((/)) = 




G = const. 



d_ 



G = const. 




G = const. 



e d tA d tB d tc d ^ 4 

The corresponding invariants are as follows: 

/i = |0'=</)2yl5C - 2{ABC + ABC + ABC), 

72 = (l)''/^uj^ABC, 

h = 2t (^ABC + ABC + ABC^ - 2C{AB + iS) - ^t^'^ucp'^ABC, 

h = B{AC-AC), 



h = A{BC - BC), 

/ \ 4 4/ • 

Iq = [ABC + ABC + ABC) + -ABC {ABC + ABC + ABC) <i?u(t)^ABC, 

V / 3 3 4 

t._l_2 

= + 2 + -//^<^^i?Cln ^- j , 

/s = + 2 + -^"'<t>ABC In ^- j , 

Ig = 2{BC - BC)A\nA - 2{AC - AC)B In B + 2{AB - AB)ClnC. 



The first invariant is the Hamiltonian of the system. 
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B. V{(P) = faV 

In this case the system psp again yields nine Noether symmetries of which six NSs 
Xi, X/j^, X^, Xj, Xs, Xg are the same as above. The additional two NSs are 



^ dt ' 3 V^dA ' "dB ' ^dCl' ^ 3' 



The corresponding invariants are as follows: 
J3 = e"* {^WABC + 12{ABC + ABC + ABC) - 8a{ABC + ABC + BAC) - Su^'^^'^ABC^ , 

ie = e""* (4a^ABC + 12{AbC + ABC + AbC) + 8a{ABC + ABC + bAC) - Sucp'^cp^ABC^ 



C. V{(t)) =aexp{l3(t)'^), kj^ -2 

The system p5|) yields five Noether symmetries of which four Xi, X4, X^, Xg are the same as 
above. An additional NSs is 

d 4(/)-'=/2 d ^ d 

''' = 'm-WT2)d^^^dc' ^ = ^°^^^t- 

The corresponding invariant is 
l2 = 2t( ABC + ABC + ABC] -2C(AB + AB] - ABC ate^^^ + , + ^ ^ 



P{k + 2) 

In this approach, we have obtained the exponential potential of the scalar field. Such exponential 
potentials are most favored in cosmology to study dark energy dynamics and fullfill many issues 
of dark energy approach, both from a theoretical point of view and in comparison with available 
observational data 23] ■ Therefore in the present model, we have obtained a solution representing 
accelerated expansion and is of immense cosmological interest. The case with A; = is interesting 
as it leads to = aexp(/30)). Such potential forms have been used a lot in phenomenological 

models of dark energy such as quintessence, phantom and quintom [24^. For this potential, we 
have acceleration solution. Further in this case, when the kinetic term has much larger pressure 
than the potential term, then the potential domination epoch is an attractor solution as long as 
the potential is flat, i.e the case /3 = 0. 
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FIG. 1: Cosmological evolution of w vs time t. The model parameters chosen as a = 1, /3 = —1, uj = \. 
Curves in various colors correspond to (red, k — 1), (blue, k = 2), (black, k — 3), (green, k — 4). 

D. Arbitrary V{4>) 

The system (llSh gives four Noether symmetries which are the same as Xi , X4 , , Xg in sub- 
section A. Therefore we don't get any new NGS and corresponding invariant in this case. 

IV. COSMOLOGICAL IMPLICATIONS 

In this section we investigate the general cosmic evolution of the model proposed in section II 
with an exponential potential given by 

y(0) = aexp(/3(/>^). (16) 

The EoS parameter w and the deceleration parameter q can be constructed analytically. Since 
w and q depend on metric functions and scalar field, their evolutionary behavior is obtained by 
numerically solving the Euler-Lagrange equations (9)-(12) for an appropriate set of the parameters 
and the initial conditions. 

The EoS parameter is constructed using the expressions of total energy density and averaged 
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FIG. 2: Variation of q vs time t. The model parameters chosen as a — 1, /3 — —1, uj = 1. Curves in various 
colors correspond to (red, k — 1), (blue, k — 2), (black, k — 3), (green, k — 4). 



pressure 



1 



p = ^io^''<j)^ + V{(l)) 



(17) 



(18) 



When A; = 0, equations (17), (18) transform to the canonical scalar field model with a rescaling of 
the field. The EoS parameter is defined as 



w 



P 
P' 



(19) 



The numerical simulations of w is drawn in figure- 1 which shows that w behaves like the phantom 
form of dark energy. This conclusion is exciting since there exists convincing astrophysical evidence 
that the observable universe is currently in the phantom phase j25|. In figures, we chose the initial 
conditions A(0) = 0.2, 5(0) = 0.1, C(0) = 0.3, </>(0) = 1, i(0) = 1, -B(O) = 1, (7(0) = 1, 0(0) = 0.3. 
Further we calculate the deceleration parameter q using the average scale factor a 



aa 
'72 



-1 



A ^ B , C_ _ (A\2 
A B C ^A' 



(4 + 



B , C^2 
B ^ C> 



(20) 



11 



Again using the numerical calculation displayed in figure-2, we show that q remains always negative 
indicating the accelerated expansion of the universe. Note that our model predicts the present value 
qo = —0.67 which is in good agreement with the astrophysical data [26 [. 



V. CONCLUSION 



In this paper, we investigated the Noether gauge symmetries of a simple extension of an old 
model proposed by Saez-Ballester in a homogenous but anisotropic Bianchi type I backgrounds. 
We solved the gauge equations and classified the models depending on potential function. One 
of the models is a generalization of the exponential families which have been used frequently in 
phenomenological models of dark energy such as quintessence, phantom and quintom. By perform- 
ing numerical simulation of cosmological parameters w and q, we demonstrated that the universe 
lies in the phantom energy dominated phase while the present value of deceleration parameter is 
compatible with the observations. 
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